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Abstract. 

In previous works, we developed a perturbative extrapolation formula to obtain 
gravitational waves at infinity from binary black hole simulations, based on the 
Regge-Wheeler-Zerilli formalism. In practice, this formula was basically derived in 
a background flat spacetime. We derive here an estimation of errors and an improved 
extrapolation formula by using the Teukolsky perturbative equation for a background 
Kerr spacetime. The improved formula requires knowledge about a mass and spin of 
the background spacetime. 
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In numerical relativity (NR) gravitational waveforms from binary black hole (BH) 
systems P |2l |3] are one of the most important output from the simulation. In P, 
we have proposed a perturbative extrapolation formula for the Weyl scalar ip^. 

’•V'M™ = rt^p(«,r) - + jM |r<^p((,r)] + 0(r-=), (1) 

decomposed in the spin-weighted spherical harmonics (_2l£m) expansion to obtain 
gravitational waves at inhnity, from a hnite extraction radius r. Here, r denotes an 
approximate areal radius, and ipf^{t,r) is the {£, m) mode of -04 at hnite radius r. We 
see that the above equation is obtained from the wave equation in a hat spacetime. It 
should be noted that gravitational waveforms h+/x are related to V’4 as '^4 = i h-^ 
and that this is true only as r —>■ cx). 

The above formula works well and has been conhrmed for example in 0 by 
comparison with a characteristic evolution code to obtain the gravitational waveform at 
null inhnity. Also, the extrapolation formula has been used in [6] for NR simulations of 
neutron star binaries. 

Here, we revisit the extrapolation formula in [1] and improve it. This paper 
is organized as follows. Section [2] gives a brief summary of the derivation of the 
extrapolation formula and the ehect on waveforms. In Section [3], the original formula 
developed in the BH perturbation approach mum is improved by use of the Weyl 
scalar -04 in NR. The improved extrapolation formula requires information of a mass 
and spin of the background Kerr spacetime. We discuss a proper choice of the mass and 
spin for the improved formula, and the diherence between the original and improved 
formulas in section 01 All throughout, we use geometric units, with G = c = 1. 


2. Perturbative extrapolation 


In the Regge-Wheeler-Zerilli (RWZ) formalism [3 [8], we consider gravitational waves. 


h. 


ihy = 




e>2, \m\<l 


-2^£m 5 


( 2 ) 


where and denote the even and odd parity wave functions, respectively. 

The wave functions satisfy the RWZ equations, 

pj2 pj2 

^ I ^ T/(e™“/odd) 

Here, r* = r -|- 2Mln[r/(2M) — 1] with a BH’s mass M, and 

are the potential and source terms, respectively. On the other hand, the NR waveforms 
are usually obtained from the NR '^4 data, '^4 = h_|_ — i hx. We should note that these 
are true only at r = cxd. 




(even/odd) 

£m 


{t.r) = S, 


(even/odd) 

im 


(t,r). (3) 
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In the analysis of the asymptotic behavior of the RWZ wave functions, we have an 
asymptotic waveform mm, 

^(^en/odd)^^^^) ^ H^^it - r*) + j dtHimit - T*) + 0(r"^), (4) 

for general (£, m) modes. Here, Hi^ is considered as the wave observed at inhnity. An 
error due to hnite extraction radii arises from the integral term on the right hand side 
of (jl]). Inverting the above relation, we have a perturbative extrapolation formula as 




(even/odd) 


Im 


= 4,(»Wo-d| I (t,+ o(r-=). 


( 8 ) 


This expression is applied to waveforms calculated in the black hole perturbation 
approach 

Next, we discuss the mode function of the Weyl scalar ip^. If the NR Weyl 
scalar satishes the Teukolsky equation [S] in the Schwarzschild spacetime, this mode 
function is described by 


r r) = Him{t - r*) + — - Ht,^{t - r*) + 0{r ^), (6) 

where the difference between and in dl]) is just a numerical factor. Here, we 
have used the Kinnersley tetrad to dehne the Weyl scalar -04 from the Weyl tensor. 
Inverting the above relation, we obtain a perturbative extrapolation formula for as 

j dtlr4'-'(t,r)] + OW^). (7) 
This is used for extrapolating gravitational waveforms in NR. 


2.1. Phase and amplitude corrections 


To see the phase and amplitude corrections by use of the above formula, we assume the 
decomposition of in (0]) as 

Hem{t - r*) = Aim exp{-iujim{t - r*)), (8) 


with the amplitude (Aim) and frequency (cvim) which are assumed to be real-valued and 
time-independent, and the phase is given as pim = c^£m(t — r*). Then, substituting ([H]) 
into (jl]), the RWZ wave functions at a hnite extraction radius are written as 

' + !)■ 


g (even/odd)/, _ 4 


1 + 


2uJimr 


exp{-iujim{t - r*)) + 0{r 


= (Aim + SAim) exp{-iuJim{t - r*)) exp{i dpim) + 0(r , 


( 9 ) 


where 6Aim and dpim are dehned as 


dAim 

-^Ira 

sin dpim 


t(l + 1) 


2oJimr 

+ + 0(r-^ 


+ 0(r-‘), 


2uimr 


( 10 ) 


The phase correction from the extrapolation formula has 0(r ^), and is the most 
signihcant correction. On the other hand, the amplitude correction will be (9(r“^) which 
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we have ignored here. When we include the 0(r“^) correction for the waveforms, the 
amplitude will be modihed. Therefore, the above analysis on the amplitude correction 
is not accurate for the coefficient. The above r-dependence is consistent with [niE], 
and has also been observed in the black hole perturbation approach [laiis]. Although 
we considered the corrections for the RWZ wave functions here, this analysis is also 
applicable to the Weyl scalar. 


3. More analysis and improved formula 


In the previous section, we have used the Kinnersley tetrad to define the Weyl scalar 
tjji which is in the Teukolsky formalism. On the other hand, the tetrad used in NR is 
different from this. Using (2.15) in [16], we discuss the tetrad dependence. 

Assuming the peeling theorem ('^4 = [r'^ 4 ]/r, "03 = = 

'ipQ = [r^'ipQ\/r^^ where the functions in the square bracket are 0{r^) for large 
r), at a finite radius we have 

1, a (7a[r^/>^^] cos^^ - 3a[rV'^^]) 


] - 


^^2 


a cos 6 [rV ’^a sin 0 + 2 M cos 6) 


+ 0(r-^), (11) 


where and are for the Kinnersley and NR tetrads, respectively. The original 
correction (see, e.g. m) of the tetrad difference in the extrapolation formula is just 

',NR1 /I 

2~V) 

Here, in order to improve the original formula, we focus only on the 0(r~^) correction. 


= 1 H-n - 


[rf^l 


( 12 ) 




M[r'ijj^^] ia[rijj^^] cos 6 


+ 


+ 0(r-^). 


(13) 


and consider the effect of the Kerr (a) term in the perturbative extrapolation. 

It should be noted that we also take care of the a dependence in the Teukolsky 
formalism. The wave function in the Teukolsky equation is - 2 ^ = (r — ia cos 
while we simply used - 2 '^ = by ignoring the a dependence in the original analysis 

from the RWZ equations im. Including the a dependence, we have 


- 24 / = (r — iacosOy 
A 


1 ^ NR _ yvk 

2 r 


iaffy^ cos9\ 

' + 0[r'-) 


+ 


r 




NR 

4 


M i a cos 9 + O(r^), 


(14) 


from the difference in the tetrads. On the other hand, the radial Teukolsky equation 
(with the Kinnersley tetrad) gives the asymptotic form of the mode function, 

Aima 


-2^1 im = Him{t - r*)r'^ + 

+ 0 (r\ (aw)^). 


- 


e(e + i) 


H imif 


( 15 ) 
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In the above calculation, we have assumed the frequency domain waveforms in flA.ip with 
a monotonic frequency, for example, Him{t — r*) = —iuHim{t — r*) in the transformation 
from the frequency to time domain analysis. This assumption works well in the case of 
quasicircular orbits because we consider the instantaneous frequency as the monotonic 
frequency. The separation constant A is given in flA.2p for aw -C 1, and we have ignored 
the 0{{auY) term. 

Comparing flTTl) with flTHl) in each {£, m) mode, 

~ / ~ 1 ) (^ + 2) Aima 

ff^mit -r ) + - --- - Hem{t -r )- - r ) 


2r 


i{i + l)r 




I'm' 


+ 0{r-\ (aw)^), (16) 


where we have considered that the NR waveform is obtained as fj^^it^r) in 


(17) 


Im 


The coefficients of mode couplings, which has a value for I' = I and = £ ± 1 

with m! = m, are calculated as 

2m 


cz = j dn_2YZ{^)cose.2YiU^) = 
cZ''"= [ dn_2YZ{n)cos9_2Ye+im{n) 


i{i +1) ’ 


'(£- 1) {i + 3){i-m + l){i + m + l) 


i + l 


{2i+l){2i + 3) 


(18) 


(see also Appendix A of [H]). It is noted that at inhnity is equal to 2H^rn{t — r*), 
and if we consider only the {£' = i, m' = m) coupling, the a term contributions cancel 
out in ffT6|) . In the following analysis, we discard the factor 1/2 in the right hand side 
of flTHl) due to the difference of the tetrads at inhnity, if necessary, this can be introduced 
again. This means that we treat 

- r*) + (^-l)(^ + 2 ) - r*) 


2 r 


i{i + l)r 


= 1 


Hul ^ w {t,r)] + 0{r 


-2 


(auf) . (19) 


I'm' 


We invert the above equation as 


H £rn{t 


= 1 


2 e'm' 


{£-!){£+ 2) ~ Aima _2 . .07 

— Himit — r ) + . .— Him{t — r ) + 0(r , (aw) ) 


= 1 - 


2r 

2M 


£{£+l)r 

l)! - Z' Y Gm kV'Sw (*. r)] 


(^-l)(^ + 2) 

2r 


I'm' 


dt[rfZL{tG)\ + £ ^ r)] + 0{r 2 , (auf) 
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-- — ^2^ + C'(r-^ (aa;)2), (20) 


where we have considered the iteration procedure which is to substitute Hirnif — G) = 
r)] for the 0{r~^) terms, and used flTSl) for the derivation of the third equality. 
Finally, we have an improved extrapolation formula 


r xfj 


im 

4 


r=oo 


(l _ Ipl (ri,Tjt,r) - + j *NS(«.r)l) 

- ^ E r)! + 0(r-^ (au,f). 

e'^e 


( 21 ) 


Here, the factor (1 — 2M/r) in the right hand side of the above equation has been 
introduced in [Ij to adjust the tetrad difference, and is consistent with the analysis 
ignoring 0(r“^) in this section. 


4. Discussion 

We have improved the perturbative extrapolation formula based on the Teukolsky 
formalism to include spinning effects in the extraction background. The improved 
formula in (|2Til is slightly more complicated than the original formula in ([71). The 
original formula will give a good result for the dominant {i = 2, m = ±2) mode of the 
waveform because there are only sub-dominant modes in the coupling due to the Kerr 
(a) term. 

On the other hand, we should be careful when we treat the sub-dominant modes. 
For example, the {i = 3, m = ±2) mode has a contribution from the dominant 
{i = 2, m = ±2) mode through the coupling. And also, we need to choose properly the 
mass (M) and Kerr (a) parameters in the improved formula of fl2T]) . Since the mass 
and spin contributions are already in a higher order with respect to 1/r, we may use 
a simple estimation, M = mi + m 2 and a = {Si -|- 5*2 -|- Lorh)/M for the aligned/anti- 
aligned cases where mi and m 2 are the mass of each BH, Si and S 2 are the spin, and 
T^orb = {mim 2 /M)yJMrorh is the (Newtonian) orbital angular momentum, respectively. 
At hrst glance, Lorb diverges for large orbital separation Torb, but the correction of the 
a/r term in fl2T|) is still small because the extraction radius should be in the wave 
zone, he., r S> For processing binary cases, we will need some complicated 

procedure to include the Kerr term. In the co-precessing frame (see, e.g., nasniEi]), 
one of the candidate of a is a = |Si -|- S 2 -l- Lorb|/Tf. More practical procedure should 
be developed for processing cases. Also, it should be noted that we have assumed the 
monotonicity of the frequency in the transformation from the frequency to time domain 
analysis. The validity of this assumption should be investigated. Recently, (l2Tll has 
been studied by setting a = 0 in [22] • 
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Appendix A. Estimation of the radiated energy, angular momentum and 
linear momentum 


To hnd the correction in the radiated energy calculation, we hrst extend the asymptotic 
waveform to the next order (r“^) in the Kerr background spacetime with mass M and 
Kerr parameter a. For simplicity, the frequency domain analysis is employed here. The 
mode function of the Teukolsky function, _2\1/ = (r — ia cos is obtained as 


-2'^emLo{r) = 


[r'’ + i I ma H— —(—3 ia + im^a + 2 Mrri) a 
2, u J V 2 


1 i {iX ma + 3 ima + 3 M) 1 A (A + 2) 

2 


r + 0(r°) 




£muj 


(A.l) 


uj 8 

where becomes the waveform at inhnity, and A is the separation constant in 

the radial and angular Teukolsky equations [9], for example, we have the {au <C 1)- 
expression [221121] as 

2 777 {-\- P -\- 4l 

A = (f + 2) (f - 1)- \ aa; + + 1) - n{^)) aV + 0{{auf ); 


+ 1 ) 

^ (f-m) (f + m) + 

{2^-l)G(2^+l) 

The energy flux from the above asymptotic waveform is obtained as 

Qauj{au — m) — X 




r = 1 + 


2uj^r‘^ 


+ Oir 


-3^ 


77100 


(A.2) 


(A.3) 


via the square of the time integration of - 2 '^imoj- Here is evaluated from the 

waveform at inhnity, H^rnuj/r. The above expression is the same as [T5| via the Sasaki- 
Nakamura equation [25]. In a similar way, we have a formula at a hnite radius for the 
radiated angular momentum [26] as 

l + 0(r-=‘))iS.„. (A.4) 




Qau{au — m) 
2uj‘^r‘^ 


The difference between and is the factor of m/oj in the frequency domain 

analysis. 

When we discuss the time evolution of linear momentum of binaries, the situation 
is slightly more complicated. Here, we focus only on the radiated linear momentum 
along the z-axis. In this case, we need to consider the mode coupling between different 
^ modes in same m modes. Formally (see, e.g. [27] for a detailed analysis in the RWZ 
formalism) 


PGr) = 






i{X — X') 24 au{au — m) — (A — A')^ — 2(A — A') 


2ujr 


+ 




poo 


(A.5) 
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where the separation constant A' is associated with i'. In contrast to the radiated energy 
and angular momentum, there is an 0{r~^) correction in the calculation of the radiated 
linear momentum. 
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